Porous Superhydrophobic Membranes: Hydrodynamic Anomaly in Oscillating Flows 
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We have fabricated and characterized a novel superhydrophobic system, a mesh-Uke porous su- 
perhydrophobic membrane with sohd area fraction $s, which can maintain intimate contact with 
outside air and water reservoirs simultaneously. Oscillatory hydrodynamic measurements on porous 
superhydrophobic membranes as a function of $s reveal surprising effects. The hydrodynamic mass 
oscillating in-phase with the membranes stays constant for 0.9 < $s < 1, but drops precipitously 
for $s < 0.9. The viscous friction shows a similar drop after a slow initial decrease proportional to 
$s. We attribute these effects to the percolation of a stable Knudsen layer of air at the interface. 



To completely describe the flow of a viscous fluid past 
a solid body, one must solve the Navier- Stokes equa- 
tions inside the fluid subject to boundary conditions on 
the solid surface |lj. These boundary conditions cannot 
be obtained from hydrodynamics, but emerge from the 
microscopic interactions of fluid particles with the sur- 
face. Consequently, they are not universal. It is well- 
established, for instance, that the commonly-assumed 
no-slip boundary condition can be violated [2 on both 
hydrophobic [3]-[6] and superhydrophobic surfaces [THTO]. 
The consequences of a shift in boundary condition from 
no-slip to partial-slip are vast. Many natural organisms 
survive simply by virtue of slip pTHT3] . Slip flows are 
expected to impact technology by enabling drag reduc- 
tion in both laminar [4] and turbulent flows [14 . This 
list goes on. 

On a conventional superhydrophobic surface [15], 
hydrophobicity combined with microscopic roughness 
causes the water surface to remain suspended above the 
solid tips, with mostly trapped air underneath [T6l [T7]. 
Since the flow is on a composite surface made up of solid 
and air, one solves the Navier-Stokes equations subject 
to no-slip on the solid elements and to slip at the water- 
air interface [2 . Thus, in a first pass analysis, viscous 
friction force on a superhydrophobic surface is found to 
be proportional to the wet solid area, ^s [18.- In this 
manuscript, we show that flow on a porous superhy- 
drophobic membrane deviates from the above picture. 
Oscillatory hydrodynamic response [19 of the membrane 
suggests that a stable Knudsen layer of gas percolates on 
the membrane, changing the boundary condition. This 
is because the porous superhydrophobic membrane struc- 
ture enables the surrounding air to move ballistically to 
the interface with little resistance — in contrast to a 
conventional superhydrophobic surface, where trapped 
gas pockets are diffusively connected to a gas reservoir 
through macroscopic distances. 

The novel system under study shown in Fig. 1 is a 
tension-dominated porous silicon nitride membrane made 



hydrophobic by silanization. The membrane has a (nom- 
inal) macroscopic area of a x a = 600 x 600 /im^ and a 
nanoscale thickness of tg = 200 nm. A matrix of identi- 
cal square pores with dimensions Ig x Ig = 10 x 10 /im^ 
are lithographically etched in the membrane. The pitch 
is Ig + /s, where Ig is the width of the solid strips in be- 
tween the pores, as shown in Fig. 1(b). This results in a 

solid area fraction $5 = 1 — t^-^-^. When a drop of wa- 
ter is placed on the porous membrane, it is supported by 
a composite surface of solid and gas (air); thus, wetting 
is not favored as shown in Fig. 1(d). 

We first characterize the intrinsic mechanical prop- 
erties of the porous membranes. In order to eliminate 
any fluidic effects, we perform these measurements un- 
der vacuum. For a tension dominated square membrane 
(a X a X tg), the frequency of the normal mode {m^n) 
in vacuum is given by ^^ = [ ^ ^^^ ^2 (^^^ + n^)]-*-^^ 
[20] . Here, as is the tension, pg is the density, and m 
and n are two integers. The in vacuo mode frequen- 
cies ^^ of a $5 = 0.34 membrane are shown in Fig. 
1(e). Here, the resonances are excited by a piezoelectric- 
shaker and detected using a Michelson interferometer 
at a pressure of 10~^ Pa. The data confirm that the 
tension-dominated membrane approximation holds well, 
even for a membrane with a very small solid fraction. 
Given that the modes are well-separated in frequency, 
each mode (m, n) can be modeled as a damped harmonic 



oscillator with effective mass Ms 



ps^stsa 



and stiffness 



^^^{m -\- n ). Relevant mechanical parameters 
for the fundamental modes of all our membranes are dis- 
played in Table 1. 

We now turn to measurements with water. The mea- 
surements are performed using a fluid cell atop the mem- 
brane as shown in Fig. 2(a). The porous membrane does 
not leak but supports intimate and continuous contact 
with both the water reservoir above and the gas reser- 
voir (ambient atmosphere) below. Using a heterodyne 
Michelson interferometer (with displacement sensitivity 
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FIG. 1: (a) Top view of a porous membrane chip [a v. a — 
600 X 600 /im^). (b), (c) Optical micrographs of $3 = 0.78 
and 0.34 membranes, respectively, (d) A drop of water placed 
on a larger membrane (a x a = 2 x 2 mm^ and $s = 0.48) 
showing the superhydrophobicity of the surface, (e) Vacuum 
mechanical resonances of a 600 x 600 \ivc? porous membrane 
with $5 = 0.34. Nearly-degenerate modes (m, n) and (n, m) 
are observed when m ^ n. Single standard deviations in the 
data are smaller than the symbols. 



TABLE I: Mechanical properties for the fundamental mode 
of the porous membranes. 
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of 100 fm/vHz around 10 kHz with 85 /iW incident on 
the photodetector), we have measured the thermal-noise 
spectra of all the membranes in their fundamental modes. 
Fig. 2(b) shows the noise spectra measured at the cen- 
ter of three membranes with different ^s- In order to 
confirm that we are working with the fundamental mode 
[(m^n) = (1, 1)], we have scanned the optical spot along 



the X and y directions and obtained mode shapes, such 
as the one shown in the inset of Fig. 2(b). Since we 
exclusively study the hydrodynamic response of the fun- 
damental mode here, we henceforth drop the superscript 
11. The top data trace in Fig. 2(c) shows all the fun- 
damental resonance frequencies in vacuum, -^^ obtained 
by driving the membranes linearly. Thermal spectra with 
water atop the membranes have provided the resonance 
frequencies ^ and linewidths ^ (bottom trace) [21]. 

We first provide a general discussion of the fluid dy- 
namics encountered in our system. We consider the out- 
of-plane (broadside) oscillations of a rigid square (a x a) 
plate immersed in a viscous fluid. We take the plate 
velocity as the real part of the complex exponential, 
Us = 5R {/7se*'^^z}, with amplitude Us- Adopting the no- 
slip boundary condition, we find the magnitude of the 
fiuidic force FfZ on the plate in the high-frequency limit 
as [231 
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and A ~ 20. The viscous boundary layer 
depends on the dynamic viscosity 



with cj) ^ I 
thickness, 5 

II f and density pf of the fluid. Mf is the so-called added 
or hydrodynamic mass, well-known from the potential 
flow theory around an accelerating solid body. Conse- 
quences of Equation (1) are as follows. Viscous energy 
dissipation is due to tangential flow on the plate, ex- 
pressed by the first term on the right-hand side. Being 
proportional to /i/, the dissipation provides a widely- used 
probe of the fluid-solid interaction. The second term on 
the right-hand side does not contribute to dissipation 
since Us^^ integrated over a cycle is zero. However, this 
term provides an independent probe of the fluid proper- 
ties (near the solid) through pf. To emphasize this, we 
write Mf = pfVf^ where Vf stands for the volume of 
fluid displaced by plate motion and depends only upon 
geometry. Indeed, it will be shown below that, in our sys- 
tem, the changes in the nature of the fluid near the solid 
boundary results in changes in both Mf and dissipation. 
Returning to the membrane oscillations, we make a 
one-dimensional harmonic oscillator approximation for 
the fundamental mode. We analyze all our experimental 
data (of Fig. 2 and Table 1) using this approximation, 
obtaining the results shown in Fig. 3. In this approxima- 
tion, the membrane has position z^, velocity Ug = igZ, 
mass Ms and stiffness Ks . We assume that Us is nearly si- 
nusoidal because all membrane resonances in water have 
quality factors Qw ^ 20 and the thermal drive has a 
white spectrum: 1^5 ^ J? {/7se*^^^}. Given that the dissi- 
pation from water dominates the overall dissipation [2T] , 
we write Msiis ^ ^e + ^w? where F^ = —KsZs is the 
elastic spring force and F^ = Ff in Eq. (1) with the ap- 
propriate parameters. Based on these considerations, we 
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FIG. 2: (a) Schematic of the measurement cell in cross-sectional and isometric views. The cell is filled with water and is 
placed on top of the membrane chip. A heterodyne Michelson interferometer probes the motion from below, (b) Thermal noise 
spectra of the fundamental mode of three different membranes with water atop. From top to bottom, $s =0.48, 0.6 and 1. The 
frequency axis is normalized with the respective resonance frequencies in water. The inset shows the shape of the fundamental 
mode for the $s = 0.6 membrane, (c) Fundamental-mode resonance frequencies in vacuum and with water atop, and linewidths 
7t(;/27r with water atop. Error bars represent the associated single standard deviations and are only shown when larger than 
the symbols. 



write a complex linear response function for the system 
as [24] G(uj) ^ [Ks - {Ms + M^)uj'^ + iaujj'K The effect 
of the fluid is embedded in two measurable parameters: 
the added water mass M,„ and the friction coefficient a 



The added water mass M^ can be determined from 
the frequency shift of the membrane mode when it is 
loaded with water. The stiffness Ks of the mode does 
not change appreciably from vacuum to water. Thus, 

^, which simplifies to Ms^v^ ^ 



mm, which results in a factor of 1/2 reduction in M^. 
Experimentally, the measured frequency increases by a 
factor of 1.2, which is close to the factor a/2 expected. For 
^s ^ 0-9, there is a significant deviation from this simple 
model: the measured M^ shows a rather fast decrease, 
eventually by a factor of 23. 

The upper, slowly increasing trace in Fig. 3(b) 
is the cycle- averaged viscous force on the membrane, 
^^^2^4'^^ Predicted using the plate model of Eq. (1) 
and accounting for the membrane mode-shape. Here, 



Ms^v = {Ms + Mu,)^, 

M^n^'^ since Ms <C M^ [27J. Figure 3(a) shows M^ as l^w is the dynamic viscosity of water; 



a function of <l>s, calculated using Ms values in Table 1 
and frequency values in Fig. 2(c). Note the two separate 
regions in Fig. 3(a) with a transition around ^s ^ 0-9. In 
order to estimate M^ from first principles, we emphasize 
that our resonator is not immersed in water [see Fig. 
2(a)] — unlike in a typical set up. There is a water 
layer of thickness t^ = 2.4 mm and density p^ atop the 
membrane, but the backside is exposed to atmosphere. 
The dominant hydrodynamic mass contribution comes 
from the the entire water layer moving in-phase with the 
membrane in the z direction [25j. This provides M^ ^ 
p^^^^f^ ^2x IQ-'^ kg [dotted segment in Fig. 3(a)]. This 
estimate is in agreement with the data of Fig. 3(a), but 
only in the region 0.9 < ^s ^ 1- Further support for our 
estimate comes when t^ is reduced to approximately 1.2 
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2ti'^VLw^ {zs'^)^ {zs'^) being the thermal amplitude of 
the membrane found from the integral of the measured 
displacement noise spectral density. It is important to 
note that, as ^s becomes smaller. Us increases. This is 
because the stiffness Ks decreases (see Table 1), the aver- 
age thermal drive force remains constant and Q^ changes 
very slowly. In the calculated plate model, the decrease in 
the wet area, a^^s^ appears to be off-set by this increase 
in Us^ thus resulting in a net increase in the drag force as 
^s decreases. The experimental cycle- averaged friction 
force is obtained from the one-dimensional damped har- 
monic oscillator model as alls ~ M^^jwUs- This force 
plotted as the lower trace in Fig. 3(b) shows a surpris- 
ing deviation from the plate model. As in added mass, 
the plate prediction agrees with the experiments when 
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FIG. 3: (a) Measured M^ as a function of $s. The line segment is the hydro dynamic mass of the entire water layer, (b) Average 
friction force and (c) the normalized friction coefficient. The plate prediction is calculated from Eq. (1) using experimental 
velocities and frequencies where needed. The normalized friction coefficient for the plate model is ~ $s- Error bars represent 
the associated single standard deviations and are only shown when larger than the symbols. 



^5^1. However, for $5 < 0.9, the drag force decreases 
rapidly, attaining a value an order of magnitude smaller 
than the plate prediction at ^^ = 0.34. 

The drag reduction on the porous membranes can be 
better assessed, if one considers the drag force per unit 
velocity: this is the friction coefficient a. Figure 3(c) 
shows the predicted and experimentally-obtained nor- 
malized friction coefficients, ^H) • ^^^ predicted value 
is proportional to ^g since the system behaves as a plate, 
but with a reduced solid area. The experimental values 
are given hy a ^ M^j^j. The data show that drag force 
for a given velocity can be reduced by a factor of 18, if one 
goes from a complete membrane (^^ = 1) to ^s = 0.34, 
i.e., a(0.34) ^ a(l)/18. 

Given that the the membranes do not leak and M^ is 
constant for 0.9 < $s < 1, we conclude that the pres- 
ence of the air reservoir does not affect the flow in this 
interval. The agrement between the predicted and mea- 
sured friction forces in the same interval provides more 
support for this conclusion. The significant deviation in 
the measured response from the plate model for $5 < 0.9 
suggests that the flow changes around $s ^ 0.9. The new 
feature of our system is its openness to air at atmospheric 
pressure. The membrane thickness, tg = 200 nm, is close 
to the mean-free-path of air, A^ ~ 60 nm. This enables 
the surrounding air to move through the membrane pores 
with little resistance. The dramatic decrease in M^ for 
^s < 0.9 can be attributed to a percolation transition: 
air bubbles localized within the well-defined pores begin 
to coalesce as ^s is decreased, eventually resulting in a 
complete gas layer, which separates the solid strips from 
the water surface. This gas layer is expected to exist in 
the Knudsen regime, with its thickness ^ smaller than its 
mean- free-path, ^ ^ A^^. 

Assuming a complete Knudsen layer at the inter- 
face, we can assess the friction reduction on a porous 



membrane with small $5. A one-dimensional model 
will suffice. We consider a large porous plate oscil- 
lating in its plane with velocity, 5R {/7se*^^^x}, under 
water with a Knudsen air layer in between the plate 
and the water. The velocity field inside the water is 
3?{/7^e-t+^^^-^-t))x} and U^ ^ Us- Since the stress 
is a continuous function of coordinate at the interface 



[z « 0), 



Here, Us^s and Uth are 



5 ' 6 

respectively the average hydrodynamic velocity and the 
thermal velocity of air molecules; pg is the density of air. 
The 1/6 factor accounts for the fraction of molecules trav- 
eling in the +z direction. Using the parameters available. 



we derive ^ 



_3_ 



The slip length [2], A 



emerges as 6 /am. at <l>s = 0.34. 

Our results might be relevant to applications. Un- 
like air bubbles on a hydrophobic surface [22 , the air 
layer in our system is stable against diffusion into the 
water because of the resistance-free influx from the air 
reservoir. Assuming that porous pipes of macroscopic 
dimensions can be manufactured, significant drag reduc- 
tion could be achieved. Several puzzling phenomena in 
bio-fluid-dynamics, including transport through and over 
bio- membranes, and propulsion over the water surface, 
may be related to the physics observed here [ISl [l3] . 
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